Spin and orbital valence bond solids in a one-dimensional spin-orbital system: 

Schwinger boson mean field theory 
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A generalized one-dimensional SU{2) x 5(7(2) spin-orbital model is studied by Schwinger boson 
mean-field theory (SBMFT) . We explore mainly the dimer phases and clarify how to capture properly 
the low temperature properties of such a system by SBMFT. The phase diagrams are exemplified. 
The three dimer phases, orbital valence bond solid (OVB) state, spin valence bond solid (SVB) state 
and spin-orbital valence bond solid (SOVB) state, are found to be favored in respectively proper 
parameter regions, and they can be characterized by the static spin and pseudospin susceptibilities 
calculated in SBMFT scheme. The result reveals that the spin-orbit coupling of SU{2) x SU{2) type 
serves as both the spin-Peierls and orbital-Peierles mechanisms that responsible for the spin-singlet 
and orbital-singlet formations respectively. 
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I. INTRODUCTION 



Looking for exotic states in transition-metal oxides is 
a fascinating problem. Spin-orbital models arose from 
the consideration of orbital degrees of freedom of d- or 
/-electron in transition metalsiSi^ When Hund's cou- 
pling, orbital anisotropy, Jahn- Teller effect, etc. are ig- 
nored, a SU{2) X SU{2) spin-orbital model is recurrently 
proposed and studied in connection with real materials, 
such as Ceo material^ spin-gap materials Na2Ti2Sb20 
and NaV2 05;^ and cubic vanadates LaVOa and YVOaii 
In this paper, we employ the Schwinger boson mean-field 
theory (SBMFT) to study the low temperature proper- 
ties of the generalized one-dimensional (ID) SU{2) x 
SU{2) spin-orbital model and reveal that the dimer 
phases in fact consist of three kinds of valence bond 
phases: orbital valence bond solid (OVB) state, spin va- 
lence bond solid (SVB) state and spin-orbital valence 
bond solid (SOVB) state. All of the three phases are 
gapped, but they can be well characterized by the spin 
and pseudospin susceptibilities, and thus are distinguish- 
able experimentally. As we will reveal, the SU{2) x SU (2) 
type of spin-orbit coupling is responsible for both the 
spin-Peierls and orbital-Peierls phenomena occurring in 
this ID model. 



This paper is organized as follows. In Section II, we 
give a detailed Schwinger boson mean-field scheme on 
the spin-orbital model. Then in Section III, we study the 
phase diagrams for two cases: (i) = 1 and T — 1/2; (ii) 
S — 1/2 and T — 1/2. The phase diagram for the former 
is novel while the one for the latter is a supplementary for 
previous works. ^'^'^ In Section IV, we deduce the spin and 
pseudospin susceptibilities. At last, a brief conclusion is 
given in Section V. 



II. MODEL HAMILTONIAN AND SCHWINGER 
BOSON MEAN-FIELD THEORY 

The one-dimensional spin-orbital Hamiltonian with the 
SU{2) X SU{2) symmetry reads 



H 



:(s 



+ y). (1) 



where S and T are two sets of spin operators which 
satisfy the SU{2) algebra with eigenvalues S and T, 
respectively, and x and y are two constant parame- 
ters. It can also be viewed as a generalized spin lad- 
der system with four-operator interactions. There are 
two special points in this model. The first is the dimer 
point D: {x,y) = {S{S + 1),T{T + 1)), whose ground 
state is perfectly dimerized and can be recognized as a 
two-fold degenerate dimerized valence bond solid.^'' An- 
other point is the FM point F: {x,y) = {-S^,-T'^), 
whose ground state possesses the maximal values for both 
^totai = I] Sm and Ttotai = J2 Tm, and is, we will see in 

the phase diagrams later, the critical point of three uni- 
form phases. The model for S" = 1/2 and T = 1/2 was 
studied both analytically and numerically, and the phase 
diagrams were proposed by several authorsi^iS*^ Except 
for the phases with conventional spin or orbital ferro- 
magnetic (FM) or antiferromagnetic (AFM) orders, the 
mixed spin-orbital valence bond state, which exhibits no 
long-range correlation and no gap near the SU{4:) sym- 
metric point S: {x,y) = (1/4, l/4)fii appears in a large 
regime in the phase diagram. The model for S = 1 and 
T = 1/2 around the point O: {x,y) = (1,1/4) was re- 
vealed to exhibit the OVB statciSiiiii 

In this paper we will employ SBMFT to study the 
phase diagram, and focus on exploring the dimer phases. 
SBMFT had been introduced to the field of strongly cor- 
related systems for a long time, and shows its merit 
in describing the spin systems at low temperatures. It 
can produce spin liquid phase and spin gap close to 
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the ground state as well as the phases with long-range 
correlationiiSiiSiiZii^ In Schwinger boson representation, 
one can introduce two boson operators a and b to repre- 
sent a spin operator S„j: 



with the local constraint 



(3) 



In this way one can introduce two types of bond opera- 
tors, 



(4a) 
(4b) 



to mimic the ferromagnetic (FM) and antiferromagnetic 
(AFM) channels of the spin S Heisenberg interactions. 
The low energy physics can be well captured by the mean 
fields, Ap = {Fmn} or — {Amn), where (...) stands for 
the thermodynamic average. Likewise, the mean fields 
for the pseudospin T are denoted by Qp and Qa for the 
FM and AFM channels, respectively. 

As revealed by Ceccatto et al.^ retaining both the 
FM and the AFM channels will lead to 



^mn^rnn (5) 



where : : denotes normal order. In the mean field scheme, 
the bond operators are decomposed as 



■FmnFmn- 



Ft 



A2 

^p: 



where the higher-order terms, (F^„ — A^) (F„i„ — Ap) 
and (Atj„ — A^) {Amn ~ A a), are ignored. We shall 
avoid applying the identity, 

'■FmnFmn - + A^^Amn ~ S , (7) 

which holds due to the constraint Eq. |j2Jl and implies 



2- Ft F 



mn-rnn--S'- for FM, 



mn^'-mn 



S-": for AFM. 



(8) 



The identity, Eq. 0, is largely violated when the con- 
straint is imposed only on average. And this is also 
true even when the Gaussian-fluctuation corrections are 
consideredi^ 

To show how Eq. JSJ works, we present a two-site 
example, H = JSi • S2, which can be solved analytically. 
After solving it by SBMFT, one can find, in the FM 
(J < 0) case, Eq. ((SJ gives ground energy /o = JS"^ 



with Ap = S and A^ = 0, and Eq. (j5J gives the same 
ground energy with Ap = S. They both are in agreement 
with the exact solution. In the AFM (J > 0) case, Eq. 
© gives fo = ~JS{S + 1) with A^ = y^S{S+l) and 
Ap^O, while Eq. ^ gives fo = -JS{S + 2) with Aa = 
\/S{S + 1). The former is just the exact energy; the 
latter approaches the exact value only in the limit S — > 
00. Thus Eq. (jHl overestimates the ground energy in 
AFM channel. We have also found that Eq. (|SJ), instead 
of Eq. can produce the exact energy of the dimer 
point in the spin-orbital model of Eq. Thus for a 

small S, the mean field scheme in Eq. lO is better than 
that in Eq. ||SJ). In practice, when solving the mean field 
equations, we found that only one of the FM and AFM 
channels can survive on a bond meanwhile for ID and 
2D unfrustrated systems. 

To perform the mean field calculation, we decompose 
Eq. into the spin S and the pseudospin T chains, 

H w '^[Js{m)Srn ■ S,„+i -t- JT{m)Tm ■ T,„+i] 

ni 

+ ^[yJT{m)+xJs{m) + Js{m)JT{'m)] (9) 

m 

where the effective couplings 



Jy(m) = (S,„ • S,„+i + x) , 
Js{m) = (T™ • T„,+i + , 



(10a) 
(10b) 



The two chains are coupled by the effective couplings, 
which will be determined self-consistently. And Eq. (jsj 
ensures a reasonable estimation of the strength of the 
effective couplings, Eq. H10|) . for both of the FM and 
AFM channels. One should notice, for 5 = T = 1/2, the 
mixed channel of S and T must be considered around the 
SU{4) point {x,y) = (1/4, 1/4),5&U-Si which is not the 
focus of this paper. In the following we apply the bond 
operators in Eq. |(SJ) for both spin S and pseudospin T 
in Eq. ©. 

There are four combinations of uniform FM and AFM 
phases: (1) S'-FM {Js = y + Q], < Q), T-FM {Jp = 
x + A], < 0); (2) S'-FM {Jp ^ y - <c)\ < 0), T-AFM 
(J5 = x -f A|, > 0); (3) 5'-AFM {Jp = y + Q], > 0), T- 
FM {Jt^x-A\< 0); (4) S-AFM {JT = y-Q\> 0), 
T-AFM {Js = X - A\ > 0). The inequalities in the 
parentheses should be fulfilled when the mean fields are 
solved. The four phases correspond to four separate re- 
gions in {x, y) plane. We omit details for these uni- 
form phases. Let us focus on more interesting cases. 
It was pointed out by Shen et al. that OVB state 
should be a preferred ground state around the point O: 
{x,y) = (1, 1/4) for S* = 1 and T = 1/2}^ This dimer- 
ization effect is also confirmed by the TMRG calcula- 
tion of orbital-orbital correlation function, and is also 
thought to be responsible for the observation on the neu- 
tron scattering experiment of YVOs.^** In fact, the gener- 
alized model of Eq. provides two other valence bond 



3 



^ o5 c>^ -o o o o 

. t , t I \ s \ I J \ 

Af Aa Aa2 
(a) (b) 

FIG. 1; Two schemes to approach the dimer phases. The 
single and double arrow lines represent T and S, respectively. 
Two parallel arrows represent FM channel, and two antipar- 
allel arrows represent AFM channel. 



solid states, the spin valence bond solid (SVB) and spin- 
orbital valence bond solid (SOVB). The three kinds of 
dimer phases have not been studied and distinguished 
within an unified framework in the previous works. To 
approach the dimer phases, we has two schemes as de- 
picted in Fig. 1(a) and (b). It turns out that the scheme 
(b) is only appropriate for exploring the SOVB phase and 
does not produce the OVB and SVB phases. So we shall 
focus on the scheme (a). 

One can find, for the scheme (a), the spin S chain and 
the pseudospin T chain are the same object mathemat- 
ically: dimerized chain with alternating FM and AFM 
bonds. For the spin S part, we can divide it into two 
sublattices A and B by assuming that the translational 
invariance is broken simultaneously. In fact, the period 
length of the chain is doubled due to the spin-Peierles 
phenomenon. We use two sets of Schwinger bosons (a, b) 
and (c, d) to represent spins on the sublattices A and B 
respectively and choose the unit cell to be composed of 
two lattice sites with one for sublattice A and another 
for sublattice B. Then the Fourier transformations read 




(11a) 
(lib) 



R,^j,6^-. (He) 
Supposed nonzero real mean fields are 



= (^2 ~ ^j^J+s)j ' (12a) 

AF = (^^[a]c,-s + b]d,^s)y (12b) 

The Lagrangian multipliers, A^'^'' and A^^"*, are imposed 
to realize the constraints on the spin S by adding the 
term in the total Hamiltonian, 



i 

+ E {^l+s^^+s + 4+sd,+s - 2S) . (13) 

3 

For a mean field approach we take A^'^'' = A^'^^ = 
Xs- After the Fourier transform, we introduce the 
Nambu spinor in the momentum space, (f>l. = 

^a[., 6_fe, cj^, d_fe, a_fc, fej^., c_fc, . One can arrive at the 

mean-field Hamiltonian for the spin S chain in compact 
form, 

k 

where 

= N [{el + y)A\ + (ei - y) Aj,] - 2NXs{2S+l), 

(15) 

and the matrix is constructed by the Kronecker prod- 
ucts of the Pauli matrices, 

k k 
+ 3a,s cos - ao®(jy®ay+ jA,s sin - ctq <^ <^ ay 

k k 
- 3f,s cos -(TQ®a^®a^+ jp^s sin - ctq® Uy® 

(16) 

with jA,s = (0| + v) Aa/2, Jks = 
{Q\~y)^F/2, jA,T = {AI+x)Qa/2, and 
iF,T — {A\ — x)Qf/2. Then we can read out the 
spectra from the poles of the bosonic Matsubara Green's 
function (8x8 matrix), 

G'^ (fc, iUn) = {ii^naz ®<Jo®<yz~- M'^) ~^ (17) 

where w„ = wk / P {n is an integer and /3 is the inverse 
temperature). Similarly we repeat the mean field calcu- 
lation on the pseudospin T chain. Both the spin and the 
pseudospin chains have two spectra, which read 

^a,. (k) = ^JXI- + + 2ujF,o.n^ (fc), (18) 

where a — {S,T), = =F, and flaik) — 

\J^a~ Ja a cos^ k. The lower quasiparticle spectrum 
'-^a,- (k) is found to be gapped at fc = ±Tr/2, which is 
related to the physical gap by the peak of the imaginary 
part of the dynamic spin susceptibility>iS. And the gap 
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is also true for half integer spin (or pseudospin) >i2iS^ By 
optimizing the total free energy 



F = eo - ^ ^ \n[nB {uJa.v {k)) {ub {l^^u (k)) + 1}] , 

k,ct.u 

(19) 

where nsiy) is the Bose-Eienstein distribution and 

eo = A^[3 (e^Al + e|A^) + 2xy + (A^ - A|,) y 
+ {Q\ - e|) X - 2\s{2S + 1) - 2\t{2T + 1)], (20) 

a set of the mean field equations are established, 



W {Ja,s,3f,s 

W [jA,T,jF,T 

X {Ja,s,3f,s 

X {jA,T,jF,T 

Y (^jA,sjF,S 

Y {jA,T,jF,T 

{ej, + y) Aa 

jA,S 

(A^ + x) 

JA,T 



25+1, 
2T+1, 
Aa, 
Qa, 

Of, 

Af 



2Ac 



JF,S 



(Ai 



= 2Xn 



JF,T 



(21a) 
(21b) 
(21c) 
(21d) 
(21e) 
(21f) 

(21g) 

(21h) 



where we have defined integrals (the symbols with 
wave decoration means dimensionless quantities, jA.a — 
jA,a/>^a,jF,a = jF,a/>^a,^a = ^a/>^a, and wc havc Sub- 
stituted the sum with integral: J2k ^ / H)' 



W{jA,ajF,a) 



X{jA,a,3F,a) 



Y{jA,a,jF,a) 



dk 

2^ 



Eh 



^jF.a 



coth 



dk (^j^ ^jF,a cos^ k \ coth 



(22a) 



2tt 



dk 

2^ 



J2 0'^'" 



, coth „ 



(22b) 



(22c) 



III. PHASE DIAGRAM 



For the transition-metal oxides, the pseudospin is usu- 
ally one half, i. e. T = 1/2, reflecting two choices 
of unfrozen Cg or t2g orbitalsii Spin-gap materials, 
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FIG. 2: (Color online) Phase diagrams for: (a) S = 1,T = 
1/2; (b) S = 1/2,T = 1/2. The uniform phases are: I for 
S-FM and T-FM; II for S-AFM and T-FM; III for S-FM and 
T-AFM; IV for S-AFM and T-AFM. The SU(4) gapless region 
is quoted from Ref.* The special points are: F - the FM point; 
D - the dimer point, S - the SU(4) symmetric point, and O - 
the representative point of OVB phase. Please see details in 
the text. 



Na2Ti2Sb20 and NaV205, are related to the S = 1/2 
case in the model of Eq. (^,^-^ while cubic vanadates, 
LaVOa and YVO3 , are related 10 5 = 1 case due to large 
Hund's coupling. ^■^•^ 

The phase diagrams are obtained by solving the mean 
field equations numerically at zero temperature. Fig. 
2(a) shows the case for 5 = 1 and T = 1/2. We 
found three dimer phases: (1) the OVB phase with 
©A %/3/2,eF = 0,Aa 7^ and Af 7^ 0; (2) the 
SVB phase with 9^ 7^ 0, Gf 7^ 0, A^ = V2 and Af = 0; 
(3) the SOVB phase with 9^ = %/3/2, 9f = 0, A^ = 
and Af — 0. As expected, the point O: {x, y) — (1, 1/4) 
lies in the OVB phase region. At the mean field level, 
the dimer phases are captured by the staggered non-zero 
AFM mean fields 9a and/or A^, which provide a per- 
fect dimer picture of spin and/or pseudospin in the whole 
phase regions. The reason for this maybe lies in two: the 
mean field treatment omits part of quantum fluctuations; 
the dimerization effect is very strong in such a system at 
least at zero temperature. In fact, TMRG has confirmed 
the orbital dimerization at the point O: {x,y) = (1, 1/4) 
since the correlation function (T„j • T„j4.i) extrapolates 
to —3/8 per bond at zero temperature, although this 
point is not exactly soluble;^ And the OVB phase is 
quite robust even when the anisotropy, Hund's coupling, 
and atomic spin-orbit interaction are taken into account 
for cubic vanadates and relevant systems. ^^'^^ We also 
guess these dimer phases may survive in two dimensions, 
although no rigorous soluble point can be referred to. A 
dimerized OVB configuration in two dimensions has been 
proposed and checked by the spin-wave theory in Ref.-— , 
which indeed exhibits lower energy than the uniform 
phases. However it is an interesting problem whether 
SOVB can survive in two dimensions. 

Fig. 2(b) shows the case for 5 = 1/2 and T = 1/2. 
In Fig. 2(b), the 5C/(4) gapless region is quoted from 
ReiA We found the SVB and OVB phases still occupy 
two areas that are sandwiched between the SOVB and 
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Xaa (q) 




FIG. 3: Static spin susceptibilities for S = 1 and T = 1/2 at 
the OVB point O: {x,y) = (1, 1/4) (solid line) and the dimer 
point D: (2,3/4) (dashed hne). 



SU{4:) gapless region. The three dimer phases with corre- 
sponding mean fields are: (1) the OVB phase with Qa = 
V3/2,eF = 0, 7^ and ^ 0; (2) the SVB phase 
with Qa ^ 0,eF ^ 0,Aa = V3/2 and Ap = 0; (3) 
the SOVB phase with Qa = V3/2, Of = 0, A^ = \/3/2 
and Af = 0. Previous works ignored the SVB and OVB 
phases. And there is a S-AFM and T-AFM phase re- 
gion from the point of view of SBMFT. 

The result here in fact suggests that the spin-orbit cou- 
pling of SU (2) X SU{2) type plays the roles of the intrinsic 
spin-Pcierls and orbital-Peierls mechanisms. The SVB 
formation in this ID model is an example that expresses 
the concept of orbitally driven Peierls stateS^ But notice 
we did not consider the lattice distortion here. For spin 
5 = 1, the SVB phase is harder to realize since the region 
shrinks. 



IV. SUSCEPTIBILITY 

The spin structures of the dimer phases are detectable 
by spin susceptibility. Let us consider the spin suscepti- 
bility (the psedospin case is handled in the same way). 
We can define the spin density waves on sublattices A 
and B (due to the rotational symmetry, we only consider 
z component of the spin), 

(9) = E (23b) 

The spin susceptibility contains contributions from both 
the intra-sublattice and inter-sublattice fluctuations. 



m,neA,B 

= XAA (q) + XAB (q) + XBA (q) + XBB [q] ■ (24) 

In Matsubara formalism, the intra-sublattice and inter- 
sublattice contributions read 

XAA {q,T) ^ {S\{q,T) S\{-q,Q)) 

-Gl,i-k,T)GfAk + q,r) 
~GlA~k,T)Gl,{k + q,T)], (25a) 
XAB {q,T) = {S%{q,T)SI,{-q,0)) 

= lf^^[2GlA-k,T)GfAk + q,r) 

-Gl,i~k,T)Gf^{k + q,T) 
-Gl-,i-k,r)GiAk + q,r)], (25b) 

respectively, where the Green's function in imaginary 
time T is related to Eq. (|17l) by 

n 

The numerical solutions of the static spin susceptibili- 
ties, XAA (g, T = 0+) and xab {q, r = 0+), for S* = 1 and 
T = 1/2 at the OVB point O: {x,y) = (1, 1/4) and the 
dimer point D: (2,3/4) are shown in Fig. 3. At the 
OVB point O, the static spin susceptibilities exhibit fi- 
nite sharp peaks at g = tt for xaa {q, t) and q = for 
XAB {q, t) reflecting the strong AFM fluctuation among 
sites in sublattice A (or B) and FM fluctuation between 
sites in sublattice A and sites in sublattice B. The finite 
peaks also indicate no long-range order exists at zero tem- 
perature. While there is no sharp peak at the dimer point 
-D, since the formation of spin singlets leads to disper- 

sionless spin spectra, ujs,- (k) — ti's,+ (k) — \J ^% ~ 3\ s- 
The spin susceptibilities at the OVB point O and the 
dimer point D characterize the whole OVB region and 
SOVB region respectively. As we go from OVB phase to 
SVB phase, the spin and pseudospin change their roles. 
In SOVB phase, the sharp peaks will disappear for both 
spin and pseudospin. Thus the three dimer phases are 
well characterized by the spin and pseudospin suscep- 
tibilities and distinguishable experimentally. A recent 
neutron scattering experiment of YVO3 has provided a 
tentative proof of the OVBiM 

V. CONCLUSION 

In conclusion, we have employed SBMFT to study the 
valence bond states as well as uniform FM and AFM 
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states in the ID SU{2) x SU{2) spin-orbital model. The 
phase diagrams are constructed by comparing the ground 
energies of the proposed possible states. In the va- 
lence bond states the translational invariance is broken in 
space, FM and AFM parameters compete with each other 
and are determined self-consistently. Specifically, two 
sublattices are introduced to approach the three dimer- 
ized SVB, OVB and SOVB phases by SBMFT. We clar- 
ified how to capture properly the low temperature prop- 
erties of such a system by SBMFT. The main result of 
the paper is that the SVB and OVB have lower energy 



than the SOVB in some regions. These consequences re- 
veal that the SU{2) x SU{2) type spin-orbit coupling can 
provide both spin-Peierls and orbital-Peierls mechanisms. 
The results still should be compared with other phases 
including those with translational invariance. Static spin 
and pseudospin susceptibilities had been calculated in 
the present theory and are available to distinguish spin- 
singlet and orbital-singlet formations in real materials. 

This work was supported by the Research Grant Coun- 
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